“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1970-09 


Application of Cohn's sensitivity theorem to 
time domain responses 


Dahlvig, Alan Lee 


Monterey, California; Naval Postgraduate School 
http://ndl.handle.net/10945/14885 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








a ee a ee ee Se eee 2 ee ee ee eee we Se 





Limited States 
Maval Postgraduate Schoo! 





APPLICATION OF COIN'S SENSITIVITY THEOREM 


TO TIME DOMAIN RESPONSES 


by 


Alan Lee Dahlvig 





pape Cmbe. gle. 





LE A AA ES A OS NE A ST A 2 rm EE EE Oe Se NP Se 





This document has been approved for public rze- 
Lease and sake; 4ts distribution 46 unLonited. 








Application of Cohn's Sensitivity Theorem 


to Time Domain Responses 
by 
Alan Lee Dahlvig 
4 


Lieutenant, United States Navy 
B.S.E.E., University of Oklahoma, 1963 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 


September 1970 


Thinner DY los 





ABSTRACT 


A sensitivity theorem by R. M. Cohn states that a linear resistive 
circuits the ratio of a fractional change in the d.c. input resistance 
to a fractional change in an internal resistance is equal to the square 
me the ratio of the current through that internal resistance to the dc, 
input current. The theorem is extended to show the sensitivity of 
input impedance to changes in internal impedances for an arbitrary network 
at all frequencies. Equations are developed which show the relation 
between Sensitivities and instantaneous power in the frequency domain. 
An extension to the time domain makes digital computer solutions 


possible. 
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L. INTRODUCTION 


In 1950, R. M. Cohn [Ref. 1 | presented a linear pana network 
sensitivity theorem which was new and significant in its approach to 
Sensitivity. Although the theorem 18 restrictive in 2ts oOricino ts some 
it is, in fact, a particular case of a more general result involving 
circuit sensitivities as shall be discussed here. 

im Chapter 2, a proof of Cohn's Theorem is presented. In Chapter 3, 
the theorem is extended to an arbitrary input impedance. In Chapter 4, 
the results of the previous chapter, which were developed in the frequency 
domain, are interpreted in terms of time domain responses. A general 
relationship between sensitivity and power relationships is developed. 
The results are then illustrated by means of an example for which a 


computer program has been written to obtain the appropriate waveforms. 





II. PROOF OF COHN'S THEOREM 


Consider a two terminal resistive network consisting of branches 


OM 2 atic lik ss. «ae b_, having indépendendent variable resistors: us ee 
n n 
The input resistance R (see Figure 1) is then a rational function of the 


ms. Consider also, that a wnit current is applied to the inpweet cua oe 


It can then be shown that R has the following property: 


2 Pe 

he eae r 
OR a OR/R nm m m 
—- = => - er ro - Ee - 
m m =m R 


where i is the current flowing in b EawatiomweGl) . which 1s eGolmes 
ig@eorem |kef. 1|, states that the ratio of the fractional change in the 
input resistance to the fractional change in an internal resistance is 
proportional to the ratio of the power dissipated in the internal resistance 
to the input power. The proof that follows is based on that of Cohn. 


Let i, denote the current flowing in, and v, the voltage across, branch 
1 d 


ree Let the nefiwork consist of q nodes designated by SUE Se iielon eles 
Strent Source connected from n, to Hae L£ the. loops are designated by 
lile,.--s), an@etive G)irection Of Current flow peumecci, Dranci sand minm@eccin 


loop assigned arbitrarily, then Kirchoff's equations may be stated as: 


is 
Se ae 
], > 5¥: = 0 Ae eee aC C25) 
edad 
j=l 
pe 
\ s 
pe i oe Lo ee ee ee (3) 
= alee J 
j=l 
DS = 1 if the branch current and loop current have the same ¢irection, 
le iiurhe: Dranchiecurrent and sloopscircrent Mave opposite directions, 
= 0 if the branch is not in the loop. 
one = | if the branch current is directed away from the node, 


II 


-l if the branch current is directed toward the node, 





O af -the branch: 1Samoet conneceted tosthe node. 


II 


fee for) j] =e 
J 

= -] for j} =q 

= 0 elsewhere. 

Define ec. = Ne to be the voltage across Be in a pOSlLtives@rrection 
and let the voltage of node n. with respect to node ae be ee = aoe Vou 
jen if t, goes from ne to n, and £, =-haet 1. ¢oes from me tonne le 

J L J J J J 1 
node q is grounded, a = Ey =R, (4) 


Since Kirchoff's laws imply the conservation of power, the power input 
must equal the sum of the powers consumed by the individual branches. 


This may be written as: 
n 

fhe dl ee = R : (5) 
ee, 


If equation (3) is differentiated with respect to te LES Lonumne 


Eimet : 
n n 
a fel 0g 
aoe es | = 5 ji geese 2 
or L, ij jd d nj) (ene or i: 1,2,...,q (6) 
m . m m 
ileal aoe 


Bet the nodes of ae be denoted by oe and ne Sicily) Ehar =e. 1S positive 
7 3 


going from ce to Sst bien 


iy 
eae i. =i hee 5 fi) ies (7) 
jj 1a jb ja j 

nN nh 

ae! ) 

Sek arte eras. = Eee : 8 
jj ab) ia 8) 


h, = 0 ae ae! |e am (9) 





and using equation (7) above, 


n nN n 

y | \ h \C E. )t ( 
es ey eee | ee eee oa) ie . 10 

ie alee) iB ie ea ja ) 

j=l fell j=l 


Since ae and oer are voltages measured with respect to node q, they 


may be written as E,‘'s. Also, a network is not necessarily a series 

i 
memnection of resistors. Thus, ae and Pee may be associated with more 
than one branch. If the sum on the right hand side of equation (8) is 


performed and all terms with the same E. are combined, 1t is seen that: 


n n q n 

i Buty 

/ Z - = hE 

L*5° 55 1 ar ee i ie L2i5°j i (11) 
sla! a al j=l 


Noe 
= 0 by the definition of 2 a : 


mest 
13 J 


This can easily be shown in a simple example (see Figure 2). 


q 
), & ~ Ey = (EY ~ E,)hy ar (Ey - Eh, aR (ED ~ E, dh, + 
k=1 
(Ey - E,)h, + (E, - EL)h, 
= Ej dh, “3 hy 2 he) ay E, Ch, - be) ae E,(-h, + h, - hy) + 
ea - hy) : (12) 


From the definition for a given earlier: 


a =a =a =a = a =e 


24 £6 


Fe oyna EV amtsren makers ar le 


and all others equal O. 
If each h, in equation (12) is multiplied by its appropriate a 
J 
such that i corresponds to the appropriate index of the Ee Tt wel ieie 


Seen that equation (11) holds. 


Comparing equation (6) and the defining equation for h,, (9), it is 


seen that Hae. wena e Tiitiswaeinoam =( 1)) 
4 








n n 5; 
Saab, =) - 3 
ee eee ar = ; (13) 
j=l j=l 
If equation (5) is now differentiated with respect to r , the result is: 
z Or i Si. 
~ ay “2 / « vac ; = 5 
OR acm |) itr,| = lee +2) se) a cs J (14) 
Or Or Te an) gal Oe ee ah) “jorene 
m m ; ; m ; m 
q=b iy jzl 
From equation (13) 
n 
Oi. 
ie 
I: 


Bsince the resistors of the circuit are independent, 


OE. { Oe em 
or by) = 


and equation (14) becomes 


OR 2 
ve = (14a) 
m 


This is Cohn's theorem for a resistive network with a unit step input 


Surrenct. 





III, EXTENSION OF COHN'S THEORE 





The expression developed in equation (14a) is ince eens but is 
limited to purely resistive circuits. Therefore, an extension of equation 
(l4a) to an arbitrary impedance is desireable. By applying a theorem of 
Tellegen [Ref. 2], it can be shown that: 


=e 


ie, Coe) 
eZ(s) _ b | (15) 
02, (s) I(s) 


The method used follows to some extent the work of Penfield [Ref. 3}. 
Tellegen's theorem states, "In a network configuration, imagine 
branch currents i such that for every node i = 0, imagine branch voltages 
v such that for every mesh Dv = 0, and for every branch let the positive 
airect ton of current be from the + to the - denoting the positive polarity 
Gt the voltage. Then Ziv = 0, a the summation is over all branches". 
Tellegen then goes on to show how for Figure (1) the theorem can be 
written as ar = a where La and v, are the terminal voltages and 
Smenments and Ly and Nh are branch currents and voltages. It is significant 
to note that Tellegen's Theorem is strictly tepological and is related in 


no way to the linearity or nonlinearity of the circuit components. For a 


(= Dr'y . -imad 
bb uv The unprime 


constant, linear system, Tellegen shows that DI 
quantities are the first state values and the primed quantities are the 
second state values of two circuits with the same topology .* 


Since the Laplace operator is a linear operator, it can be applied 


pod tellegen's theorem can be written as: 


*It is unfortunate that Penfield chose to use the term "state values" 
Since it can easily be confused with the "state variables" often used in 
Sarcuit analysis. As used here, “state values'' refer to the values of the 
voltages and currents of the branches of two different circuits. These two 
circuits have the same topology but not necessariiy the same circuit 
elements. 
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= Ne 
fi, MSV SS) = J, (Sov Gs) 


(16) 
t b 
Or 
)i(s)v, (s) = /V Cs) 1 Cs) Cl) 
t b 
Subtracting (17) from (16) yields: 
7 oF 
2 (2, CV's) - Thsyv.Gs)] =) [ve@rt,@) - v1) (18) 
C b 


Unless otherwise stated, all calculation in the remainder of this 
paper will be in the Laplace domain (s - domain). 
Consider circuits as shown in Figure 3. 
Hem sci cule Ca) 


Pi@icwee ) 1GCll td ble 


Vis) I(s) V(s) 


I(s) 
Pont Ey L(s) a) TGs eae nC) 
areuit V.¢s) I (s) Vs) “+ 6V, (s) I (s) ++ ore) 


Using equation (18), 
oe SS) = | Geet TAGS) Ue = Y 11, Ly, (s) “+ 6V, (s) | ~ 
v,(s)EL, (s) + 61, (s) I} 
ea) 1<S) =) {v, (s)61, (s) . 1, (86%, (9)} 

b 


V¢s) = 4 (s)L(s) and for small changes 


6V (s) = 2 (s)61, (s) te I (s)64 (s) 


(20) 

Multiplying (20) by I (s) and rearranging will yield: 
1, (s)8V, (s) - V, (8)61, (8) = Te(s)S2 (8), (21) 

From the two states in Figure 3, it can be seen that: 
bye = I(s)Z(s) 2) 

and Be Sai ks eons) IbZiG=) 02 (5) | 


I 


Eis Zcs)e tr Ligi62(s) +2 (s)oi(sjot 6165 )64(5 ees 


Substituting equation (22) into (23) shows that: 


Meayoccs et Ziteyolts) + 51(s)o2(s) — 0 (24) 
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Tao -E(s)61(s) - I(s)61(s)62(s) 


= BG ES tees) 62 Gy ae (25) 
Meri(s)6Z(s) << E (8) lien OAs) k= 2. (5) ete. 
ACN OaCE) = -E,(8)81(s) (26) 
pmpstitutinge (26) and (21) into (19), it is found that: 
Tae) >) I, (s)62, (s) (27) 
b 
| . ae _V 32s 
emace Z(s) is a function of Zz (s); CZs) =) ee 6Z, (s) 
oa 465 


fitting this result into (27) results in 


Po) aS 6Z, (s) = 2 oe) 
b 





I (3), 
a Lier Sree 
ae Gay Se L ee “ey 
b 
b 
32 (s) -[2 aa a 
oZ, (s) Cs) ; 
This is Cohn's theorem in the s - domain for any impedance. 


Bes 15° Set sequal Co ju, the results can bewebtaincdsinwener equcne, 


domain: 


ri » (JX) Z 


Zao ~| (29) 


Zz, (iw) rata 


Mie right hand side of (29) may be written as a transfer function as a 


magnitude and a phase angle. 


IT (jw) , 

eee Sa | 48 Ww) 

wee ooo: (30) 

Ee eee fai) 2 2 

62 (ju) _| _b oe 

8Z, (ju) LiCje) J” T’(w)/20@) (31) 
If Z Cj) =R, 

ae = = T° (w)/20@) (32) 


vz 





Tf 2 (jw) = job, 


one 2 
Saale = T (w) /20 Ww) 


and 
of hie) = jwt* (w) /28 (w) = w2* (w) /28 (w) eae (33) 


ia Z (jw) = 1/jwc, 


and 


7 Z 
Sttay = C1/ja)T" @) 2e@) =? (2e@) = W/2 (34) 


Thus, the magnitude of the sensitivity changes as the product of taney 
2 
or T (w)/w and its phase is advanced or retarded by 1/2 when the element 


of interest is an inductor or a capacitor rather than a resistor. 


i 





IV. COHN'S THEOREM IN THE TIME DOMAIN 


Equation (28) expresses Cohn's theorem in the s ~- domain for any 


input current I(s) applied to an impedance Z(s). 


eZ(s)_ _ ples 
3Z, (s) L(s) | 





Meus equation 1S easy to use if the input current I, and the interna! 


element current I, are known as functions of s. The input current wave 


b 
shape will be known for most application and its transform to the s - domain 
will usually not be difficult. However, for any but the most elementary 


Seecircuits, the internal element current I (s) can be found only through 

a series of tedious calculations. In contrast, there are several good 
computer programs which will calculate i (t). It is, therefore, desireable 
to find a direct means of determining the sensitivity as a function of 

time rather than by taking the inverse transform of (28). For the class 


of input currents given by (39), as derived later, it can be shown that: 


2m 

0Z(s) a “a k-2n_2 

82, (s) (eae 2 Bee : (35) 
k=0 


where n and m are integers defined by (39). This leads to the following 


mecurtcs in the time domain: 


2n 
-2n+ 
22 Sg Ba ee C00) (36) 
Xp k k- 2n+j b 2b 
k=0 
where the symbol * is used to indicate convolution, and 


i} meet ee x iS an haeawe tor (L) 


= O if X, is a resistor (R) 
= -] ifx, is Veapacitor s(t G) 


A very interesting simplification occurs when ee = 6(t), a unit impulse 
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Serttiat (Ss) = A(se 


OE(s = 
TROVE) na (37) 
Se, /z, = e, Cedai, Ce) (38) 
b 


Equation (37) shows that the sensitivity is related to the instantaneous 
power in the frequency domain. 

In attempting to develop the time domain equivalent of equation (35), 
it would appear to a caSual observer that a simple application of the 
Principal of convolution would apply. That is, X(s)Y¥(s) = £[x(t)sy(t) ] 
fgeethus, X(s)/Y(s) = LLx(t) xl/y(t) |. However, this is not true since 
ie o(t) 1 Z oy For example: 


Let Ys) 


II 
” 


VCC otic oun len 


Now let 1/Y(s) ips 
IVSrtramstorms sO mwinitestepsin tne 
time domain and a unit step is 
not the reciprocal of the unit 
doublet. 

It would appear uSeful to try to determine if there is a direct way 
Hemobtain the inverse transform of 1/Y(s) to get a time function which can 
be convolved with x(t). There does not appear to be any such simple 
Operation. 


pou alternate possibility is to limit the consideration to input 


current wave shapes which can be manipulated easily so as to yield a direct 


solution. Consider an input current of the following form: 
os" N(s) 
Gs) = = = ie) (39) 
% j 
[yore 
j=0 
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g = any real number 
n = any integer +, -, or 0 
b. = all real numbers (all but one can be Zero). 


This allows the input to include impulses, steps, ramps, sinewaves, decaying 


Svea: | Sinh at 


or increasing exponentials, : ; ees : 


Considering the input current specified in equation (39), equation (28) 








becomes 
nm , 
Wale) D(s)I Cs) 
TSS 2 ce) 
b N (s) 
2m 
ean X lone 
=p s }, bs ICs) 
k=0 
2m 
ee 2 a) 2 
= b I 
LA Nees ge 
k=0 
2m 
_ ae 2 
) (a, s I (s) ; (41) 
k=0 
faigeeresult is equation (35). 
Consider the impedance of the elements. If the element is an 
moaumetor, then 2, = sL and equation (28) may be written as 
O4(s) _ Ts) 
Te, ca ) 
KANG) ) = (b—) (42) 
OL I(s) 
b 
If the element is a capacitor, then Z. = L/sC, and equation (28) becomes 
Il 
az(s)__ 1 (pS? (43) 





a(1/C,) ~ 8 ues 


“The subscripts are changed because the b's are now the coefficients 
of the squared polynomial. 
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If the element is a resistor, equation (28) is 


I (s), 2 
OAS) yn ae | . 
SR ere) 3 (44) 


meus in general the sensitivity can be written as 





d2(s) _ _j Tene 
=e (45) 
Ox I(s) 
b 
minere $j = FLaif Xx is an inductor 
= O if x 1S 2 esse 
= -] if a is a compliance 
Combining (45) and (41) yields the result 
2m 
o4(s) _ 23 y¢ keen 2 
OX) : are JD 
k=0 
2m 
Sa k-2nt+j\_2 
=D (aye) (46) 
k=0 
im -can be Shown that multiplication by s in the s - domain is equi- 
walent to differentiation with respect to t in the time domain, that 
fivasion by s in the s - domain is equivalent-to integration in the time 
domain, and that multiplication of F(s) in the s - domain is equivalent to 


convolution in the time domain.® Ties, DEO Ene sexpencnt Ol eo ingamterm oF 
equation (41) is greater than zero, the individual terms of (46) become 


k-2n 


ae 5 Kean oe 


If the exponent of s is negative, the term of (46) becomes 


ay ae Li, (t) ‘k i, (t) ] 


“See, for instance, Holbrook (Ref. 4) pp. 66, 68 and 108. 


Ly 





OZ 
The number of integrations is equal to k-2n. The sensitivity, (2 ; 
b 


then becomes a sum of integrations and/or differentiations of the 
@onvolution of the current through the element of interest azten 1tselt, 


and equation (46) can be written as: 


2m 
Rel ecg) 
ee ) a ea Ge) il (t) | 
Ox -k k-2n+j b i ob 
b ou 
=0 
ims is equation (36). 


To show the power relations of equation (37), Ohm's law is applied 





(see Figure 4). Again, unless otherwise stated, all computations are made 
in the s - domain. 
je Gay) er 
a2(s)_ | p 9 
32) (s) LT (s)J 


Ete > t (s)Z2(5) 
& 


OE(s) = Hoes) 2 ey 
o2(s) = TS) 
& 
Thus : BE(s)/I(s) 1, (s)42 
az, (s) . LZ (s) 
sets) _ 16) 





& 
AE (s ; I (s) Is) ; E (s)L, (s) as 
82, (s)/Z, (s) 7 Zz, (s) re 7 LS) 


If the input current is limited to an impulse in the time domain, Le) 


will be unity and equation (42) becomes 


QE(s 7 
aA (21S) | Eo ty 6 om 
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Tf a circuit is excited with an impulse current (see Figure 5), the 
ratio of the voltage across any element Co Ehe inpubecurrent may be 
Wratten aS a transfer function such that: 

E, (s 
<> 


L,(s) 





=e 1S) , (50) 


MEethe circuit is now excited by a current source equal to the current 
through the element of interest, the new element voltage, E(s), will be 

E(s) = T(s)T, Cs) ‘ Gib) 
Memens1s true since the transfer function T(s) is a function of the network 
maagemot Of the input current or voltage. The original input current was 
an impulse so that equation (50) can be written as E(s) = T(s). 
emestatuting this into (31) will yield 

E(s) = E (LG) : (52) 
The input current can be written in the time domain as i (t), and E, 6s) £, (s) 
can be written as e, (t) te i, (t). Thus, (52) becomes in the time domain 

Oe) Me ce, (t) * i, (te) : (53) 
Equation (53) is the right hand side of equation (38). Thus, in order to 
measure the sensitivity of equation (38), the following procedure can be 
used. Excite the Oricind) circuit’ by an impulse ancmea louie i, (t). 
Mine excite tle crcuit with i, (t). The voltage across ae element is then 
equal to the sensitivity. 


There are several interesting variations of equation (28). For 


aa | 


M@etance: since Y(s) =Z (s) 
ais) 1 aus) nf fo 
32, (s) 27 (s) 02 (s) 22 @ I(s) | T(s)z(s) 4 
aY(s) _ ee ae (55) 
32, (s) obs a) 


1g 











oY(s) a its oe prean BA 
az, (s)/Z, (8) ~ LE(s) J LEG) J OW) 
feet = EY is substituted imto equation (28) 
I (s)-2 oo. (s)Y. (sa 
y) SS 





oZ(s a _ 
6Z, (s) a | eee ) a 


y (aern) : Poe 
v- (s)2%, (s) LE(s) 3 





- 2 
Remembering that Z = Y a the partial derivative yields 6Z =-6Y/Y . 


Ptecmne this into (56) results in 





oh Al) Se “oy 














ay (s)  - EG) ae 
For another variation, Z = a is used. Thus 
i ee le OMG) ee ee 
oY, (s) ae oY, (s Tee) Es) 
az(s) . [bt (58) 
oY, (s) L i(s) 1 
7 (5 a, Y (s)E, ts) _ EB (s)1(s) a 
ok (s)Y (s} a. ie 


feet) is restricted to a delta function, then equation (GS)% (56), and 


(57) may be written as 


Ch) ? 

3z, (Ss) I (s) (60) 
oY(s) _ 

az, (s)/Z(s) ~ BE NE Soe 

eS ee | 

Sy, (s) = EB. (s) (62) 


If i(t) is restricted to a delta, equation (58) and (59) may be written as 
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Of (Sse 
ay, (s) Hes) (63) 


=o 04/05) eee , 
oY, (s)/¥, (s) IL (SE Gs) (64) 


If the denominator of the right hand side of these equations are 
restricted to the forms of equation (39), time domain functions similar to 
Baeat of equation (47) can be obtained. For any sensitivity calculation 
however, the impulse deans eeanecion gives the simplest and most direct 


mesult. 
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Figure 1 


Iwo Terminal Network with Varied le Ty 


Used for Development of Cohn's Theorem 


ee 








eee 2 


Five Branen Network sedl ac 


Demonstrate Cohn's Theorem. 
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\)) a sie meee ene | 





Db) State. 


Figure 3 


Iwo Circuits with the Same Topology 


But with Different Internal Impedances. 
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Fipure 4 


Two Terminal Networks Used for Developing 


the Power Ke latuonss 


Efva constant inpul Cumnentwis mequo means. 
current source is connected to the input 
terminals. If a constant voltage source is 
required, a voltage source is connected. 


In ertner case b= ae 
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Figure 5 


Two Terminal Network Used for Developing 


the Lime Domain Sclution. 
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APPENDIX A 


Example Using a Parallel RLC Circuit 


For the circuit in »figume 6.) let Ro — 1,03y = Ges eee 
Calculate the sensitivity of the input impedance with respect to each of the 
Circuit components for a delta input current. Also, determine the sensitivity 
of the input voltage with respect to 6Z/Z for each component with the same 
delta input. 

An impulse was approximated by a triangular pulse 8.0(10°’) seconds wide 
and PS (10°) amps at its center. dhis input current was used as data for 
an ECAP solution to arrive at the currents through and voltages across the 
Senponents. 

The data obtained from the ECAP solution was used as input data to 
Subroutine SEN which is in the Computer Program section of this paper. Since 
outputs were obtained from ECAP every 50 calculations, the time increment 
for SEN was eke The final time was a). The delta was chosen 
to give it unit area so COEF is 1.0. The Laplace transform of a delta is 
a constant. Thus, SN and IPR were set at zero. The IELEM and ITYPE values 
were determined by the type of element and type of sensitivity required. 


The results are shown in Figures7 through 12, 


Zz 








Figure 6 


Parallel REC Cirewmt Used for Sample.Calculations: 
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